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Diffusion-limited aggregation (DLA) assumes that particles perform pure random 
walk at a finite temperature and aggregate when they come close enough and stick 
together. Although it is well known that DLA in two dimensions results in a ramified 
fractal structure, how the particle shape influences the formed morphology is still un- 
clear. In this work, we perform the off-lattice two-dimensional DLA simulations with 
different particle shapes of triangle, quadrangle, pentagon, hexagon, and octagon, re- 
spectively, and compared with the results for circular particles. Our results indicate 
that different particle shapes only change the local structure, but have no effects on 
the global structure of the formed fractal cluster. The local compactness decreases 
as the number of polygon edges increases. 
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I. INTRODUCTION 



Various aggregation mechanisms^ have been proposed to theoretically investigate disor- 
dered growth under non-equilibrium conditions, such as nanoparticle and colloidal aggregation]^ 
among which the diffusion-limited aggregation (DLA) has been intensively studied by many 
researchersP^ The original DLA model proposed by Witten and Sander^ performs a ran- 
dom simulation at a finite temperature on a two-dimensional square lattice. In this model 
a particle is initially fixed at the origin, and more particles are then released one by one 
and perform random walk in the space until they become close enough and stick on the 
central cluster. The above DLA procedure generates a statistical self-similar structure whose 
scale-invariant properties can be described by the fractal geometryP^^I Power law scaling of 
the two-point correlation function was discovered in the initial work by Witten and Sander^ 
and the fractal dimension of a cluster formed on a two-dimensional lattice is 5/3 regardless 
of the lattice geometry.^ 

DLA is a successful abstract model for qualitatively understanding irreversible aggrega- 
tion of ramified fractal structures observed in many experiments.^^ However, there still 
exist many other experimental observations which beyond the explanation given by DLA. 
For instance, aggregation can result in regular dendrite fractal structures, such as mag- 
netic a-Fe203,^ fractal assembly of copper nanoparticles,-- and the snowflake structure in 
nature P which may have direct connections with the experimental observations^^ that 
particle shape plays an important role in the formation of those structures. Consequently, 
some simulations^"^ have studied the effect of particle anisotropy on DLA morphology. 
Liu et al .^studied the influence of the monomer anisotropy to the DLA structure and they 
concluded that anisotropic monomers still lead to fractal patterns. Mohraz et al^ inves- 
tigated colloidal rod aggregation in three dimensions by both experiments and simulations 
and found that the fractal dimension increases with increasing rod aspect ratio. Menshutin 
and Shchur^Sl found that different degrees of monomer anisotropy result in clusters with 
different fractal dimensions and the noise-reduction level can change the morphology of the 
clusters. Nevertheless, no studies have been done to investigate the influence of particle 
shape to the formed morphology of DLA. 

In this work, we perform the two-dimensional off-lattice DLA simulations with different 
particle shapes of triangle, quadrangle, pentagon, hexagon, and octagon, respectively, and 
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analyze the local and global properties of the finally formed fractal structures, and compare 
with the results for circular particles. Our results indicate that different polygon particles 
lead to different local structures, but they have negligible effect on the global behavior of 
the fractal structure. In addition, the local compactness decreases as the number of polygon 
edges increases. The paper is organized as follows: the simulation and analysis methods are 
described in Section [II], the results are shown in Section III, followed by conclusions given 
in Section HVl 



II. METHODS 

In this section, we describe our DLA simulation method, the global structure analysis 
methods, the skeleton algorithm identifying the main branches of fractal structure, and the 
orientational order parameter calculation characterizing the local structure. 



A. Simulation method 

Our two-dimensional DLA simulations with circular particles are the same as the original 
one by Witten and Sander^ except that we perform off-lattice instead of on-lattice simula- 
tions. Additional simulations were conducted with one of the five polygons as the particle 
shape: triangle, quadrangle, pentagon, hexagon, and octagon. The radius of the circle is a 
and all the polygons are regular ones whose circumscribed circle has the same radius of a. 

In all simulations, a seed particle is initially fixed at the origin. At each step, a free particle 
is released at a random position with a distance d = df + do from the seed particle, where 
df is the distance between the seed particle and the outer-most particle on the cluster and 
d is the distance between the free particle and the outer-most particle on the cluster. The 
released free particle is allowed to translate in any directions with a random displacement 
generated from a uniform distribution [—dt, dt] and to rotate with a random angle generated 
from a uniform distribution [— 5 T ,8 T ]. 

After each movement, the distance d n between the free particle and the nearest particle 
on the cluster is calculated. The free particle stops moving only if d n is smaller than a 
critical distance d c , and in the circular case, its position is adjusted so that its center has 
a distance of 2a from the center of the nearest particle. For polygons, the orientation is 
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also adjusted along with the distance so that its edge overlaps with the closest edge of the 
nearest particle. After the free particle sticks on the cluster, a new free particle is released 
and the above procedure repeats. 

The distance do between an outer-most particle on the cluster and the initial position 
of a released particle is related to the particle concentration of a real system. The cutoff 
d c in the simulation corresponds to the effective range of the adhesive interaction between 
particles. The amplitudes of d t and 5 r reflect the system temperature. The periodic boundary 
condition^ was applied to avoid the escape of particles from the simulation space, and the 
neighbor list algorithm 26 was adopted to accelerate the simulations. In this work, we set 
a — 1, d c — 3, dt — 0.5, 5 r = 0.5 in radian, the simulation box to be a square with the side 
length L = 1000, and do = 100 to make sure that the DLA process is in a low concentration 
condition. Nine independent runs have been performed and each run contains M = 10000 
particles. 

B. Fractal dimension 

Historically, the fractal dimension of DLA was initially evaluated by Witten and Sander^ 
through the radial distribution function (RDF) g(r) of the cluster. According to the sta- 
tistical self-similar property of a DLA cluster, the RDF has the form g(r) oc r Df ~ D where 
Df is the fractal dimension of the cluster and D is the dimension of the Euclidean space in 
which the cluster is embedded. Later on, most experiments and simulations have utilized a 
more convenient relation between the number of particles and the distance from the origin 
to determine the fractal dimension Df. 27 The number N(l) of particles inside a circle with 
radius / and centered at the origin can be written a d 28 * 29 * 

N(l) = k x (l/a) D t, (1) 

where the prefactor ko, related to the lacunarity of the cluster,^ is different for clusters 
formed by different aggregation mechanisms.^ Oh and SorenserP*' argued that the relation 
is valid only when I is larger than a critical length l c ~ 10a. According to Eq. [TJ a fractal 
cluster with a larger Df has more particles within the same radius from the origin, so the 
fractal dimension D[ can be used to characterize the global compactness of a cluster. The 
prefactor ko can be used to characterize the local compactness of a cluster, as illustrated by 
our later local structure analysis. 
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C. Skeleton algorithm 



The skeleton of a cluster can be computationally identified to study its large scale 
properties,^ and from the skeleton of a cluster we can easily obtain the number of main 
branches rib- The main branches are those branches whose length, which is the number of 
particles from the tip of the branch to the origin, is larger than a critical length L c , which 
is comparable to the cluster size. 

In a cluster, there exist many branches whose length is comparable to the cluster size, but 
not all those branches are intrinsically different because some of them share lots of particles 
and only differ in a few particles. If the number of shared particles exceeds an critical 
value S c , all these branches are regarded as in the same main-branch class and only one of 
them is chosen to represent this class. According to the algorithm developed by Schwarzer 
et aZ.^for the DLA cluster growth process, there exists a parent-child relation between two 
neighbor particles, in which the particle joining the cluster later is called "child" and the 
particle joining the cluster earlier is called "parent". By applying this relation repeatedly 
we can identify all the branches from a tip particle to the seed particle at the origin. After 
all branches of the cluster are identified, we can obtain the main branches of the cluster as 
follows: 1) the branches whose length is larger than a critical branch length L c are picked 
up from all the tip branches; 2) the similarity S, defined as the number of particles shared 
by two branches, are calculated for each pair of branches; 3) the branches with S larger 
than a critical similarity S c can be considered as belonging to the same main-branch class 
and only one of them is chosen randomly to represent this class. The skeletons obtained 
with the above algorithm for the DLA clusters with different polygon particles are drawn 
in Fig. [T] We can see that the selected main branches are not necessarily the longest in a 
main-branch class since they were randomly chosen. 

With the main branches of the cluster determined, we then calculate the angle 9 between 
two neighboring main branches and the direction of each main branch is determined by fitting 
the main branch with a straight line. The distribution of angle 9 quantifies the rotational 
symmetry of the skeleton. The direction of the fitted straight line is mainly determined by 
the particles shared by all the branches in the same class, so the random selection of the 
main branches does not influence noticeably the calculated rotational symmetry. 
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FIG. 1. DLA clusters with different polygon particles of a) triangle, b) quadrangle, c) pentagon, 
d) hexagon, e) octagon, and /) circle. The red lines depict the calculated skeletons. 
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D. Orientational order parameter 



In contrast to the fractal dimension defined to quantify the global structure, the orien- 
tational order parameter is used to characterize the local symmetry of a DLA cluster. The 
nth orientational order parameter $ n is defined as^l 

1 - JVt . 1 . Zi 

^ = ^E-E em ^> ( 2 ) 

* »=1 % 3=1 

where N t is the total number of particles in the cluster, Z{ is the number of nearest neighbors 
of particle i, and 6y is the angle of the vector from particle % to particle j with respect to 
a fixed vector do, which was chosen in our calculation to be the unit vector parallel to the 
y axis. According to the definition, $ n equals to 1 if the local structure has a perfect n-fold 
symmetry. To observe the change of the orientational symmetry with respect to the length 
scale, we extend the definition of the orientational order parameter to be distance-dependent: 

1=1 1 j=\ 

Eq. [3] only differs from Eq. [2] in the number of particles in the second sum z' i: which is 
now over the particles in the range < r, rather than all the nearest particles. With the 
extended orientational order parameter, we can analyze the cluster geometry at different 
length scales. 

III. RESULTS AND DISCUSSION 

A. Fractal dimension 

In this section, we present our calculation of the fractal dimensions for the DLA clusters 
formed by different polygon particles through the relation between N(l) and / defined in Eq. 
[IJ The particle numbers N(l) versus the normalized radius from the origin l/a are shown in 
a log-log plot in Fig. [2] with only the linear part ranging from 25a to 148a. All the lines in 
Fig. [2] are parallel to each other and the interceptions of these lines decrease with increasing 
number of particle edges. According to Eq. [TJ the slopes of these lines correspond to the 
fractal dimensions of the clusters and the interceptions correspond to the prefactors ko- The 
fitted fractal dimensions and their standard deviations are listed in Tab. IB We can see from 
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FIG. 2. Log-log plot of the number of particles N as a function of the normalized distance l/a 
from the origin. 



both Fig. [2] and Tab. [I] that the clusters formed by different polygon particles have similar 
fractal dimensions but different prefactors. 

In Tab. |TJ the DLA clusters with triangle, quadrangle, and hexagon particles have the 
same fractal dimension of about 1.69, consistent with the fractal dimension of 5/3 for the 
on-lattice simulations performed on a triangle, square, or honeycomb lattice.^ Furthermore, 
our results show that other particle shapes of pentagon, octagon, and circle, which have no 
corresponding lattice structures, also have similar fractal dimensions. This result indicates 
that the shape of particles has no noticeable influences on the global structure of a two- 
dimensional DLA cluster, but influences significantly the prefactor ko defined in Eq. [T], 
which characterizes the local compactness of the cluster. The prefactor kg will be studied 
in detail in subsection IIIID1 
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TABLE I. Fractal dimensions Df and their standard deviations a. 

Particle D{ a 



triangle 1.67 0.04 

quadrangle 1.69 0.07 

pentagon 1.69 0.07 

hexagon 1.69 0.04 

octagon 1.75 0.08 

circle 1.73 0.07 



B. Cluster skeleton 

As described by Schwarzer et al.f^ the number of main branches nb obtained by the 
skeleton algorithm is also an important property characterizing the global structure of the 
cluster. They have also shown that, for the DLA cluster with circular particles in two dimen- 
sions, the number of main branches takes a constant value of nb = 7.5 ± 1.5, independent of 
the cluster size. In this work, we obtained the average values of nb for triangle, quadrangle, 
pentagon, hexagon, octagon, and circle, respectively, as listed in Tab. [TTJ which are all close 
to the value 7.5 ±1.5 reported in Ref. 31. 

TABLE II. Number of main branches rib and their standard deviations a. 

Particle nb o 



triangle 6.89 1.05 
quadrangle 6.22 1.09 
pentagon 6.89 1.05 
hexagon 6.75 1.28 
octagon 7.25 0.71 
circle 6.78 0.97 



Schwarzer et a/.^has shown that, in a two-dimensional DLA, the increasing rate of free 
space is the same as the increasing rate of space screened by the dangling branches aside 
the main branches, so the number of main branches is a constant during two-dimensional 
aggregation. Consistently, our calculated skeletons shown in Fig. [T] have similar structures 
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FIG. 3. Distributions of angle 9 between two neighboring main branches in the skeletons of the 
DLA clusters with various particle shapes. 



for different polygon particles. The distribution of angle 9 between two neighboring main 
branches was then calculated to characterize the rotational symmetry of the skeleton. Fig. [3] 
indicates that the distributions of angle 9 are similar for all the clusters formed by various 



polygon particles, and the average angle values are listed in Tab.|III[ All these results indicate 
that the particle shape has no effect on the skeleton of DLA cluster, neither on the number 
nor on the structure of the main branches. 



C. Orientational order parameter 

The local structure of the DLA cluster is important for studying the effect of different 
physical conditions on the growth process. Mandelbrot!^ studied in detail the lacunarity 
distribution of cluster, which is related to the compactness of the cluster, and found that the 
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TABLE III. Angle 6 between two neighboring main branches and their standard deviations a. 

Particle 9 a 
triangle 0.91 0.40 
quadrangle 1.01 0.46 
pentagon 0.91 0.44 
hexagon 0.93 0.43 
octagon 0.87 0.37 
circle 0.93 0.40 



lacunarity distribution is different at different length scales. In our work, the orientational 
order parameter was calculated to quantify the local structure of the DLA cluster. We 
calculated with Eq. [3] six orientational order parameters^ $ 3; $ 4; $ 5; $ 6) $ g; an d $10 as a 
function of distance for all simulations. The order parameters $ n as a function of distance 
r are shown in Fig. |4| 

As can be seen from Fig. |4j the value at the smallest r of each curve for the nth orien- 
tational order parameter is always 1 for a polygon with the n-fold symmetry because the 
nearest neighbors of a particle are always arranged with the orientation determined by the 
polygon shape. All figures in Fig. [4] indicate that, at the length scale a < I < 5a, the 
orientational symmetry is apparently different for the clusters formed by different polygons, 
but at the length scale I > 10a, all clusters have the same isotropic orientational symmetry. 
The local structure of a cluster is determined by the competition between the randomness 
due to diffusion and the orientational symmetry due to polygon shape around a particle. 
The influence of polygon shape on a cluster's structure decreases rapidly with increasing 
distance. As shown in Fig. [4j when I < 5a, the particle alignment is mainly determined 
by the particle geometry; when a DLA cluster is growing, the number of new growth sites 
which can accept new particles increases exponentially and distribute isotropically around 
the particle at the origin, so when the length scale / > 10a, the randomness due to diffusion 
dominates and the cluster structure is independent of the particle geometry. 
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FIG. 4. Orientational order parameters with different orders versus radius r for the clusters formed 
by particles with various shapes. (a)Three-fold, (b)four-fold, (c)five-fold, (d)six-fold, (e)eight-fold, 
and (f) ten-fold. 
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D. Local compactness and prefactor 



The above results for orientational order parameters indicate that the particle shape only 
affects the local structure of the cluster, and in Fig. [2] it is shown that the prefactors ko 
are different for different particle shapes. Therefore, there should exists a direct relation 
between the local structure and the prefactor ko- According to the original definition of 
fractal structure,^ the fractal dimension Df is calculated by measuring the volume of fractal 
structures embedded in a D-dimensional Euclidean space. The volume of fractal structures 
is obtained by counting the number M(r) of balls with radius r needed to cover the fractal 
structure. For an ideal fractal structure, the self-similarity is satisfied in all the length 
scale, 1 so the numbers of counting balls with radius rx and r 2 , respectively, should satisfy 
the following relation 

M(r 2 ) = M(r x ) x ( ri /r 2 ) Di , r x > r 2 . (4) 

Form Eq. [4] we can define the average number of balls with radius r 2 in a ball with radius 

r% as 

M(r 2 ; ri ) = ^4 = (rx/ra)^, r x > r 2 . (5) 
The above relation can be easily extended recursively to the case with three radii 

M(r 3 ; n) = M(r 2 ; n) x M(r 3 ; r 2 ), n > r 2 > r 3 . (6) 

The quantity M(a; I) is equivalent to N(l) defined in Eq. [TJ A detailed study by Oh and 
SorenserP^ concludes that the self-similar relation does not satisfy at all length scales, rather 
there exists a critical length scale of about 10a, when I > 10a, the DLA cluster has different 
self-similar property for different aggregation kinetics, but no differences were found for 
/ < 10a. Recently, Heinson et alP^ studied the inertia tensor of many different fractal 
clusters and pointed out that the prefactor k is related to the cluster's morphology. In 
addition, our results for the orientational order parameter indicate that the particle shape 
only affects the cluster structure at finite length scales. These results suggest that there 
should exist a cutoff distance l c for the recursive relation in Eq. [6j and the relation between 
the number of particles and the distance from the origin can be written as 

Nil) = N(l c ) x M(l c - 1) = N(l c ) x (l/Q D *, (7) 
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where N(l) and N(l c ) are the numbers of particles within radii I and l c , respectively, as 
defined in Eq. [TJ Compared with Eq. [TJ we obtain the relation between k and N(l c ) as 

k Q = N(l c ) x a D </l?. (8) 

By fitting the simulation results, we found that Eq. M is roughly satisfied for all cases when 



l c = 5a. The data of ln(JV(i c = 5a) x a D{ /1® S ) and ln(fc ) are compared in Tab. IV The 
prefactor ko for each case was determined by the linear fitting of N(l) vs I averaged over 
nine independent simulations. It is clear that ko corresponds to the local compactness of 
the cluster when I < 5a and the local compactness decreases as the number of polygon 
edges increases. These results are in agreement with Heinson et a/.'s conclusion^ that ko 
characterizes the anisotropy of the cluster which only exhibits in local structures. 

TABLE IV. Prefactor values for the cluster with various particle shapes. 



Particle ln(N(l c = 5a) x a Df //f f ) ln(fe ) 


triangle 


0.312 


0.464 


quadrangle 


-0.233 


-0.114 


pentagon 


-0.326 


-0.335 


hexagon 


-0.465 


-0.398 


octagon 


-0.738 


-0.769 


circle 


-0.79 


-0.82 



IV. CONCLUSION 

In conclusion, we have done a series of two-dimensional off-lattice DLA simulations with 
particle shapes of different polygons, and the global and local structures of the clusters have 
been studied in detail. By analyzing the fractal dimension, we conclude that the geometry 
of the particles has no effects on the global structure of cluster. We have also shown that the 
DLA clusters formed by different polygon particles have the same number of main branches 
and the same skeleton symmetry. The results of the extended orientational order parameters 
show that the particle shape influences the local structure of the clusters at a finite length 
scale, but the effects decay quickly at larger length scales. The prefactor connecting the 
number of particles and the radius from the origin quantifies the local compactness of the 
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DLA clusters, and our results for the prefactor indicate that the local compactness decreases 
as the number of polygon edges increases. Therefore, we conclude that the particle shape 
only affects the local structure of a two-dimensional DLA cluster, but has no effects on its 
global structure. 
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